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We present isovector nucleon observables: the axial, tensor, and scalar charges and the Dirac 
radius. Using the BMW clover-improved Wilson action and pion masses as low as 149 MeV, we 

£f) achieve good control over chiral extrapolation to the physical point. Our analysis is done using 

three different source-sink separations in order to identify excited-state effects, and we make use 

f^**) of the summation method to reduce their size. 
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1. Introduction 

Lattice QCD calculations of nucleon structure observables are now entering the stage where 
they are being seriously confronted with experiment, although resources are not sufficient for full 
control over all systematic errors. Thus, when predicting unobserved quantities, it is useful to know 
which sources of systematic error must be more carefully controlled in order to achieve agreement 
for experimentally observed quantities, and which sources of error are already well under control 
when using standard techniques. 

We report on two experimentally measured observables, the isovector Dirac radius {r\) v and 
the axial charge g&; and two predictions, the tensor and scalar charges, gj and g$. Proton matrix 
elements of the vector current are parameterized by two form factors, 

(p(P')\qf q \p(P)) = u(P>) (fF?(Q 2 ) + ^F«(Q 2 )) u(P), (1.1) 

where A = P' — P and Q 2 = —A 2 . The isovector Dirac radius is defined from the slope of the 
isovector = F" — Ff at zero momentum transfer: 

F?(Q 2 )=F?(O)(l-l(r 2 yQ 2 + 0(Q 4 )). (1.2) 

In terms of experimental observables, (r 2 ) v is related to the difference between proton and neutron 
charge radii, the former of which has a 7a discrepancy between measurements from electron- 
proton scattering [1] and a recent result using a precise measurement of the Lamb shift in muonic 
hydrogen [2]. Lattice QCD calculations with a precision of a few percent could help to resolve this. 

The axial, tensor, and scalar charges all have similar definitions from neutron-to-proton transi- 
tion matrix elements at zero momentum transfer: 

{p\uY^Y 5 d\n) = gAUp'fysUn, {p\uo^ v d\n) = gjUpO^ u n , {p\ud\n) = g s u p u n . (1.3) 

The axial charge is a key benchmark observable, since it is a naturally isovector quantity (and 
thus not requiring disconnected diagrams to calculate it), measured via forward matrix elements, 
and it is also well-measured experimentally via beta decay of polarized neutrons. 

The tensor and scalar charges have not been measured experimentally, but it has recently been 
shown [3] that they control the leading contributions to neutron beta decay from new (beyond the 
Standard Model) physics, and thus they provide a useful input to the analysis of experimental data. 

2. Methodology 

The main results presented are from calculations performed on ten Lattice QCD ensembles 
using 2 + 1 flavors of tree-level clover-improved Wilson fermions coupled to double HEX-smeared 
gauge fields [4]. We also compare with results from earlier 2+1 flavor calculations [5, 6, 7] 
using four ensembles with unitary domain wall fermions [8, 9], as well as five ensembles using a 
mixed-action scheme with domain wall valence quarks and Asqtad staggered sea quarks [10]. The 
ranges of parameters used in these calculations are summarized in Tab. 1 and Fig. 1 . 

On every ensemble with Wilson fermions, we compute nucleon three-point functions using three 
different source-sink separations T . When matrix elements are computed using the traditional ratio- 
plateau method, the asymptotically dominant excited-state contaminations arise from transitions 
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Action 


a (fm) 


m % (MeV) 


L x /a 


Lt/a 


T/a 


# meas 


Wilson 


0.09 


317(2) 


32 


64 


10, 13, 16 


824 


Wilson 


0.116 


149-356 


24, 32, 48 


24, 48, 96 


8, 10, 12 


762-10032 


Domain wall 


0.084 


297-403 


32 


64 


12 


4216-7056 


Domain wall 


0.114 


329(5) 


24 


64 


9 


3192 


Mixed 


0.124 


293-597 


20, 28 


64 


9 


2176-5024 



Table 1: Lattice actions and ranges of parameters used: lattice spacing, pion mass, spatial and temporal box 
size, source-sink separation, and number of measurements. 




Figure 1: Spatial box length L x and pion mass 
for the full set of ensembles. The areas of 
the circles are proportional to the number of 



0.1 0.2 0.3 0.4 0.5 0.6 

mj (GeV) measurements made on each ensemble. 

between the ground state and the lowest excited state, and these decay as e _A£T / 2 . The summation 
method requires combining more than one source-sink separation, but yields improved asymptotic 
behavior, with the leading contaminants to forward matrix elements decaying as Te _A£r [11, 12]. 

This set of ensembles allows for control over, or study of, various sources of systematic error. 
In roughly decreasing order of the level of control that we can achieve: 

Quark masses For all ensembles, the strange quark mass is near the physical value. Our smallest 
m u d corresponds to a pion mass of 149 MeV, which is just 10% above the physical pion mass. 
This allows for either a direct comparison of the m % = 149 MeV ensemble with experiment, 
or a comparison after a mild chiral extrapolation to the physical pion mass. 

Excited states Using the ratio-plateau method with three different source-sink separations allows 
for clear identification of observables where excited-state contamination is a problem. These 
three source-sink separations can also be combined using the summation method to get another 
result that may be less affected by excited states. 

Finite volume In general, finite- volume effects are expected to be small with m n L > 4. Furthermore, 
we can perform controlled comparisons between the 24 3 x 48 and 32 x 48 Wilson ensembles 
near m n = 250 MeV, where the spatial volume is changed while other parameters are fixed. 

Finite temperature At our smallest pion masses, the lattice time extent is shorter than the typically 
used L t = 2L X , and these ensembles may be particularly susceptible to thermal effects. On 
the other hand, the three different time extents L t /a € {24,48,96} used for Wilson ensembles 
near m n = 250 MeV are useful for identifying possible problems. 

Discretization The use of different lattice actions and different lattice spacings allows for a consis- 
tency check, but this set of ensembles is insufficient for taking a continuum limit. 
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Figure 2: Isovector Dirac radius (r^) v , as determined from dipole fits to Fi(Q 2 ). The two experimental 
points both use the PDG [1] value for (r|)", and (r|) p is taken from either the PDG or from the result 
from measurement of the Lamb shift in muonic hydrogen [2], Left: Results from ensembles using the three 
different lattice actions. Wilson action points are taken from the middle source-sink separation. Right: The 
full set of Wilson action results. Measurements using the three source-sink separations and the summation 
method are slightly displaced horizontally. The points corresponding to the smallest source-sink separation 
are placed at the measured value of m\, except for the ensembles with m n w 250 MeV and m K w 350 MeV, 
where the different volumes are displaced horizontally and the dotted lines indicate the approximate measured 
values of m|. 



3. Results 

For each of the observables, two plots are shown. The first includes the data from all three 
actions. In order to show results using similar techniques, the middle source-sink separation on the 
Wilson action enembles is shown for this comparison. The second plot shows the dependence on 
source-sink separation and the summation result for each of the ten Wilson ensembles. 

Figure 2 shows the isovector Dirac radius. This is extracted from measurements of the isovector 
Dirac form factor F^(Q 2 ) via dipole fits F{'(Q 2 ) ~ jjT^m 1 ^' ^ e ^ rst P^ ot sriows a broad con- 
sistency among the three lattice actions in their region of overlapping pion masses, although the 
unitary domain wall data are systematically a bit higher than the mixed action data. There is a 
gentle rise as m % approaches the physical value, but all of these data significantly undershoot the 
experimental results. The second plot shows that excited-state effects are responsible for a large 
part of this discrepancy with experiment. There is a general trend for the data on each ensemble to 
increase with source-sink separation, and on the lightest ensemble the summation point is near the 
experimental points. Chiral extrapolation of the summation data to the physical pion mass yields a 
result consistent with experiment [13]. As shown in Ref. [13], calculations of the isovector Pauli 
radius, magnetic moment, and quark momentum fraction behave similarly: results for the 149 MeV 
ensemble monotonically approach the experimental value with increasing source-sink separation, 
the summation point agrees with experiment, and a chiral extrapolation of the summation data to 
the physical pion mass yields results in agreement with experiment and having a smaller statistical 
uncertainty than the 149 MeV summation point. Further study will be needed to obtain full control 
over excited-state effects, but with presently available data we consider the summation method as 
our best approach for reducing their size, and it succeeds in producing agreement with experiment 
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Figure 3: Axial charge g^. The experimental value is from Ref. [1]. See caption of Fig. 2. 
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Figure 4: Tensor charge gj- The physical pion mass is indicated by the vertical line. See caption of Fig. 2. 



for these four observables. Comparing the summation data on the four different space-time volumes 
at m n ~ 250 MeV, we see that there is no statistically significant dependence on spatial volume or 
time extent. The summation point on the fine ensemble at m % = 317 MeV is also consistent with 
those on the nearby coarse ensembles, indicating the absence of discretization effects at this level of 
precision. 

The axial charge gA is shown in Fig. 3. Again there is broad agreement among the different 
actions. There is no clear dependence on m n , and (using the middle source-sink separation for the 
Wilson ensembles) the data undershoot experiment by about 10%. Looking at the second plot, we see 
that, on the lightest two ensembles, increasing the source-sink separation moves the data away from 
experiment. The opposite behavior is seen in a subset of the four m % « 250 MeV ensembles. For the 
three with L t /a> 48, increasing the source-sink separation moves the data toward experiment and 
the summation points are consistent with experiment. In contrast, the fourth with L t /a = 24 behaves 
similarly to the lightest two ensembles, albeit with larger statistical uncertainty. This suggests that 
the decrease of gA with source-sink separation is caused by the influence of thermal pion states, 
since the three ensembles that show this behavior all have small m n L t . 

Figures 4 and 5 show the tensor charge gj and the scalar charge g$, respectively. Measurements 
of the latter are much noisier; note the significantly larger range of the vertical axis in Fig. 5 
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Figure 5: Scalar charge g$. The physical pion mass is indicated by the vertical line. See caption of Fig. 2. 



compared with Figs. 3 and 4. Neither g$ nor gj shows a clear statistically significant dependence on 
source-sink separation, so using the middle source-sink separation should give better results than 
was the case for (r\ ) v and g&- At m n Pt 250 MeV, there isn't a significant dependence on the spatial 
volume, but the ensemble with small L t shows a dependence on source-sink separation for both 
gs and gj. This behavior isn't seen on the ensembles with the smallest pion masses, so it is likely 
that they aren't strongly affected by thermal states. From the middle source-sink separation on the 
m n = 149 MeV ensemble, we get g$ = 1.01(27) and gr = 1.049(23). Extrapolation to the physical 
pion mass using either the coarse Wilson ensembles, or a global fit to all ensembles, yields results 
consistent with the lightest ensemble [14]. 

4. Conclusions 

The importance of near-physical quark masses for nucelon structure calculations is illustrated by 
the isovector Dirac radius, where the rise toward experiment is only seen at our lightest pion masses, 
and for the axial charge, where new behavior was seen only below m n Pt 250 MeV. In addition, it is 
essential that excited-state effects can be identified, as shown clearly for the isovector Dirac radius. 
This shows up even more dramatically for the isovector quark momentum fraction [15, 13]. 

More study of excited-state effects is required, and this may require large computing resources, since 
as the source-sink separation T is increased to reduce excited-state contamination, the signal-to-noise 
ratio is expected to decay as e -( m N-~2 m ") T [16]. We have identified finite-temperature effects as a 
possible source of the discrepancy with experiment for the axial charge. 

As we obtain a better understanding of systematic errors, predictions of nucleon properties 
using Lattice QCD become more credible. Calculations of the nucleon scalar and tensor charge will 
provide useful input to searches for new physics. 
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